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Abstract 


The  buckling  loads  of  8-ply  graphite-epoxy  delaminated 
cylindrical  panels  were  determined  experimentally.  The 
analysis  included  two  different  ply  orientations,  two 
different  aspect  ratios  (length/chord),  two  different 
delamination  sizes,  and  one  set  of  boundary  conditions; 
clamped  along  the  top  and  bottom  edges  and  simply 
supported  along  the  vertical  sides.  The  experimental  test 
results  were  compared  to  the  linear  bifurcation  loads 
obtained  from  the  STAGSC-1  finite  element  computer  code 
and  the  buckling  loads  of  panels  with  square  cutouts  which 
were  also  obtained  from  STAGSC-1. 

The  experimental  testing  was  accomplished  by  the  Air 
Force  Flight  Dynamics  Laboratory.  Experimental 
difficulties  were  encountered  due  to  the  inability  of  the 
test  apparatus  to  evenly  distribute  the  load  through  the 
simply  supported  boundary  conditions.  The  uneven  load 
distribution  resulted  in  nonsymmetric  buckling  patterns, 
however,  the  buckling  loads  showed  good  correlation 
between  analytical  and  experimental  values.  All 
experimental  bucKling  loads  were  less  than  the  analytical 
value  predicted  by  STAGSC-1.  The  amount  of  decrease  in 
load  depended  on  ply  orientation,  aspect  ratio,  and 


•lei  am  mat  ion  size. 


Equations  (9)  and  (10)  are  expanded  for  an  N-layered 


laminate  as  : 
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where  zk  and  zk_i_  are  defined  in  Fig.  4. 

By  substituting  Equation  (8)  into  Equations  (11)  and 
(12)  and  recalling  the  stiffness  matrix,  Qj_j  is  constant 
within  a  lamina,  the  integration  is  simplified  by  pulling 
the  stiffness  matrix  outside  the  integration.  The 
stiffness  matrix  is  summed  over  each  layer  and  Equations 
(11)  and  (12)  are  rewritten  as: 
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With  these  assumptions,  the  following  strain-curvature 


relation  can  be  written  for  a  laminate: 
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where  the  superscript  'O'  indicates  the  mid-surface 
strains,  z  is  the  distance  from  the  panel  mid-surface,  and 
K  are  the  mid-surface  curvatures.  Substituting  Equation 
(7)  into  Equation  (4),  the  stress  in  the  kth  layer  can  be 
expressed  in  terms  of  the  laminate  middle  surface  strains 
and  curvatures: 
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Once  the  stresses  in  a  layer  are  known,  the  resultant 
forces  and  moments  are  obtained  by  integrating  the 
stresses  in  each  lamina  through  the  total  laminate 
thickness.  The  force,  Nx ,  is  expressed  as  a  force  per 
unit  length  and  moment,  Mx ,  expressed  as  a  moment  per 
unit  length,  are  calculated  from: 
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Ql6  =  (Qii~Qi2“2Q66  )  sine  cos3e  +  ( Q 1 2“Q22  +  3<266  *  s  ln  3#  cos# 


Q26  =  (Qii-Qi2_2Q66 > sin3s cose  + (Qi2-Q22+Q66> sln^cos3e  <5-Cont- > 

Q*6 6  =  (Qii+Q22-2Qi2-2Q66)sin2ecos2e+Q66  (sm4e+cos40) 

The  stress-strain  relation  defined  in  Equation  (4) 
represents  the  relation  for  the  k^  layer  in  a  laminate: 


To  use  Equation  (6)  and  expand  the  stress-strain 
relations  to  a  multilayered  laminate  (Fig.  4)  certain 
assumptions  are  made.  First,  the  intralaminar  bonds  are 
considered  perfect  so  that  ro  lamina  can  slip  relative  to 
another.  Second,  the  Kirchof f-Love  hypothesis  is  valid 
and  requires  that  sections  normal  to  the  panel's  mid¬ 
surface  remain  plain  and  normal  to  the  mid-surface  after 
bending  [ 16 ] . 


where  E]_  and  E2  are  Young's  moduli  in  the  1  and  2 
directions  respectively.  Poissons'  ratio,  Uj_j  is  the 
ratio  of  transverse  strain  in  the  j  direction  when  loaded 
in  the  i  direction.  G^2  is  the  shear  modulus  in  the  1-2 
plane.  The  values  are  further  defined  in  the  following 
relations : 

y12  =  '  E2/El 

V21  =  ~  E 1 / E 2  (3] 

G1 2  =  T12/712 

The  preceding  relations  defined  the  strains  and 
stresses  in  only  the  principal  material  directions. 

Common  use  of  composites  have  the  lamina  axis  oriented  at 
some  angle,  6,  to  the  structural  axis.  In  this  case  the 
stresses  are: 
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where  Qij /  the  transformed  reduced  stiffness,  is 
defined  as: 

Qll  =  Q^cos4#  +2  (Q12  +  2<266  )  sin2  9  cos 2? +Q22sin4? 
Ql2  =  ( Q 1 1  +Q22-4Q66  )  5 in^ cos  ^0  "'"Q12  (sin^^+cros^©) 
G22  =  Qnsin46 +2  (Q^2  +  2G66  *  si-n2®  cos2?  +Q22  cos4? 


10 


cannot  maintain  equilibrium..  The  lower  path  represents 
the  occurrence  of  buckling  and  the  usual  linear 
load-stress-strain-deflection  relations  do  not  hold. 
However,  the  lower  path  represents  equilibrium  in  the 
structure . 

Classical  Lamination  Theory 

In  order  to  understand  the  behavior  of  a  composite 
laminate,  it  is  necessary  to  understand  the  basic 
constitutive  relationships.  References  [16]  and  [17] 
present  an  in-depth  development  of  these  relations. 

Before  a  complete  laminate  made  up  of  N  layers  can  be 
analyzed,  the  relations  for  a  single  lamina  are 
presented.  The  stress-strain  relations  for  the  individual 
laminae  are  similar  to  those  of  an  orthotropic  material  in 
plane  stress.  This  relation  in  principal  material 
coordinates  (Fig.  3)  is  given  as: 


where  Q^j  j_s  the  reduced  stiffness  matrix  defined  in 
terms  of  the  engineering  constants: 

Qll  =  El/i_  y12  U21 

Ql2  =  V 12  E2/1-Ui2u21  =  u21El^~  u12  ^21  <2> 


Q22  =  E2/1-u12u21 
066  =  g12 
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End-Shortening 


Load-Deflection  Curve 
Figure  2 
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Bifurcation  Buckling 

During  the  past  25  years,  the  study  of  the  stability 
of  shells  has  become  one  of  the  most  investigated 
structural  problems  [14].  As  a  consequence  of  these 
investigations,  many  theories  of  buckling  are  available, 
but  none  has  yet  become  universally  accepted.  Yet,  the 
results  of  these  studies  have  contributed  to  the 
understanding  of  many  aspects  of  shell  stability. 

Buchert  [15]  defines  buckling  as  "collapse  of  a 
structure  under  stress  or  loads  that  are  less  than  those 
causing  material  failures."  With  this  definition,  the 
term  buckling  is  used  synonymously  with  the  bifurcation 
point  throughout  this  thesis.  The  concept  of  buckling, 
and  the  bifurcation  point,  can  best  be  explained  with  a 
plot  of  load  versus  end-shortening  of  the  panel  [14]. 

Fig.  2  is  a  typical  plot  of  load  versus  end-shortening  for 
a  thin  she  1 1 . 

For  small  loads  and  deflections,  the  curve  is  a 
straight  line  representing  the  equilibrium  state  of  the 
shell.  As  the  load  is  increased,  the  bifurcation  point  is 
reached.  At  the  bifurcation  point,  there  are  two 
theoretically  possible  paths.  The  upper  path  implies 
buckling  is  prevented  and  the  usual  load-stress-strain- 
deflection  relations  hold  true.  However,  this  state 
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xtytz  =»  structural  coordinates 
u,VfW  =  displacements 
t  =  thickness 
R  =  radius  of  curvature 
C  **  chord 
L  »  length 
©  «■  ply  orientation 


Panel  Notation 


Figure  i 


[0/90] 2s-  Three  different  size  de laminations  were 
included;  0  inch  (no  delamination),  2  inch  circular,  and 
4  inch  circular.  The  panels  were  manufactured  with  8-ply 
graphite-epoxy  with  all  de laminations  placed  at  the 
midsurface  and  located  at  the  geometric  center  of  the 
panel.  A  total  of  24  panels  were  tested  which  represents 
a  duplication  of  each  aspect  ratio,  ply  orientation,  and 
delamination  size. 

The  material  properties  were  determined  experimentally 
for  a  given  ply; 


El 

=  20. 

.97  X 

106 

psi 

e2 

=  1. 

.53  X 

106 

psi 

g12 

= 

.815  X 

106 

ps 

V12  =  0.300 

The  panel  boundary  conditions  incorporated  in  the 
experimental  test  device  provided  a  clamped  top  edge 
(u  =  v  =  w  =  w;X  =  0) ,  clamped  bottom  edge 
(u  =  free,  v  =  w  =  wfX  =  0);  and  simply  supported  on 
the  vertical  sides.  See  Fig.  1  for  a  description  of 
panel  dimensions  and  displacements. 
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work  of  Yin  and  Fei  [13]  refined  the  analysis  to  include  a 
circular  delamination  in  a  thin,  circular  plate.  Their 
results  indicate  a  panel  under  pure  compression  will  only 
buckle,  while  a  panel  loaded  under  both  compression  and 
bending  will  experience  buckling,  snap-through  (separation 
of  the  delaminated  layer  from  the  main  body  of  the  plate) , 
and  catastrophic  failure  if  the  load  is  sufficient  to 
attain  the  critical  value  for  growth  of  the  delamination. 


Purpose 

The  purpose  of  this  thesis  is  to  experimentally 
determine  the  buckling  loads  of  cylindrical  composite 
panels  with  circular  delaminations  positioned  in  such  a 
fashion  that  the  shell  skin  and  delamination  buckle 
together.  Limiting  values  on  the  experimental  loads  are 
determined  from  the  linear  bifurcation  load  calculated  by 
the  STAGSC-1  Finite  Element  Computer  Code.  Three  major 
variables  are  included  in  the  study:  ply  orientation, 
aspect  ratio,  and  delamination  size.  The  results  are  used 
to  further  the  understanding  of  composite  buckling 
behavior  and  increase  the  data  base  available  to  users  of 
composite  materials. 


Scope 


The  experimental  panels  consisted  of  two  different 
aspect  ratios,  1(12  X  12)  and  1  1/2  (12  X  8),  and  two 
different  ply  orientations,  [0/-45/ +45/90] s  and 
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compared  to  theoretical  values.  These  investigations 
continued  [4,  5]  as  both  experimental  and  analytical 
techniques  were  refined,  resulting  in  a  broader,  more 
reliable  data  base  that  provides  comprehensive  set  of 
buckling  load  data  for  a  wide  range  of  panel  aspect  ratios 
(length/chord)  and  ply  orientations.  To  continue  the 
understanding  of  the  buckling  problem,  researchers  have 
begun  examining  the  effects  of  different  variables  on  the 
buckling  problem.  Boundary  conditions  are  one  variable 
that  has  a  substantial  effect  on  the  buckling  load  of  the 
panel  [6,  7].  These  investigations  indicated  that  the 
buckling  load  increases  as  degrees  of  freedom  are 
restrained  at  the  boundary  condition,  resulting  in  a 
stiffer  panel.  Another  variable  considered  in  recent 
studies  is  panel  imperfections  [8].  These  imperfections 
range  from  eccentricities  caused  by  the  manuf actur ing 
process  or  test  fixture  installation  to  cutouts  that 
remove  a  percentage  of  the  panel's  area  [9]. 

Recently  numerous  articles  have  dealt  with  the  effects 
of  a  delamination  on  an  axially  loaded  laminated  plate. 
Analytical  approaches  to  the  problem  [10,  11,  12] 
developed  techniques  to  determine  a  flat  panel's  axial 
load  capacity  for  a  variety  of  boundary  conditions, 
delamination  lengths,  and  plate  and  delamination  thick¬ 
nesses.  However,  these  techniques  were  limited  to  a 
one-dimensional,  across-the-width ,  delamination .  The 
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INSTABILITY  OF  COMPOSITE  CYLINDRICAL  PANELS 

I .  Introduction 

Background 

The  benefit  of  composite  material's  high  strength  to 

weight  ratio  make  it  an  excellent  material  for  use  in 

modern  aircraft  structures.  The  USAF 1 s  commitment  to  the 

use  of  composite  aircraft  structural  parts  has  resulted  in 

the  Air  Force  Aeronautical  Systems  Division  recently 

awarding  a  contract  to  McDonnell  Douglas  Corp.  for  the 

development  of  a  computerized  system  to  control  the 

autoclave  curing  process  [1].  Also,  NASA  is  currently 

engaged  in  developing  and  demonstrating  the  efficient  use 

of  composites  in  airframe  structures  on  large  transport 

aircraft  [2].  This  widespread  use  of  composite  materials 

has  accounted  for  an  increase  in  the  amount  of  analytical 

and  experimental  research  on  the  behavior  of  composites. 

However,  one  area  of  research  still  in  its  early  stages  is 

the  effect  of  a  delamination  or  debonding  on  the  overall 

stability  of  a  laminated  panel. 

Initial  development  of  procedures  to  experimentally 

determine  the  buckling  load  of  curved  composite  panels  was 

done  by  Wilkens  [31,  who  also  presented  a  percent 

knockdown  factor,  ( ( 1-exper imenta 1  buckling  load)  X  100 j, 

theoretical  buckling  load 

to  show  the  decrease  in  experimental  buckling  loads  as 


Further  simplifications  to  the  integrations  are  possible 
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by  recalling  the  strains  and  curvatures  in  Equations  (13) 
and  (14)  are  mid-surface  values  and  not  functions  of  z. 

Hence,  they  too  can  be  removed  form  the  integration  and 
the  moment  and  force  equations  become: 
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where  A^j  are  the  extensional  stiffnesses,  B j_ j  are  the 
coupling  stiffnesses,  and  D^j  are  the  bending  stiffnesses. 
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Delamination  Panel  Buckling/Postbuckling  States 

There  are  three  distinct  states  to  the  buckling/ 
postbuckling  problem  of  panels  containing  de laminations 
defined  by  Fei  and  Yin  [13]:  bifurcation,  snap-through , 
and  catastrophic  failure.  The  bifurcation  phenomena 
results  from  the  panel  being  subjected  to  pure  compression 
and  the  delaminated  layer  maintains  full  contact  with  the 
main  body  of  the  plate.  With  the  introduction  of  bending 
and  the  proper  ratio  of  bending  moments  to  the  compressive 
load,  snap-through  of  the  delaminated  region  v.  L 1 1  occur. 
Snap-through  is  defined  as  separation  of  the  delaminated 
layer  from  the  main  body  of  the  plate.  A  continued 
increase  in  the  load  past  the  snap-through  states  results 
in  an  increasing  energy  release  rate.  When  the  energy 
release  rate  reaches  the  critical  value  for  growth  of  the 
de lamination ,  the  delamination  will  grow  and  catastrophic 
failure  will  occur. 
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In  their  analysis  of  the  postbuckling  problem,  Fei  and 
Yin  [13]  developed  three  non-dimensional ized  expressions 
for  bending  strain  versus  membrane  strain,  delamination 
radius  and  energy  release  rate.  The  expressions  considered 
a  homogeneous,  isotropic,  circular  plate  with  a  centrally 
located  circular  delamination  subjected  to  axis  symmetric 
edge  compression  and  bending  moment  along  its  boundary. 

The  three  relationships  are: 

p  =  vh/2 t 

-  bending  strain  versus  membrane  strain  (19) 

a  =  [12(1  -v2)t  ]  1/2  a/h 

-  non-dimensionalized  delamination  radius  (20) 

Ga4  =  |l2(l-u2)  |2  G*a4/Eh5 

-  non-dimensionalized  energy-release  rate  (21) 

where : 

M  -  curvature  to  the  plate  middle  surface 

h  -  depth,  measured  from  top  surface,  of  delamination 

f  -  compressive  membrane  strain  in  delaminated  layer 

v  -  Poisson's  Ratio 

a  -  delamination  radius 

G*  -  material  fracture  toughness 

E  -  Young's  Modulus 


Their  results,  presented  in  table  and  graph  form, 
relate  the  three  non-dimensionalized  terms  to  determine 
which  buckling/postbuckling  state  the  plate  will 
experience . 

This  thesis  is  concerned  with  determining  the  initial 
buckling  load  of  a  panel  with  a  delamination  located  at 
the  midsurface.  The  results  of  Fei  and  Yin  are  used  to 
confirm  that  due  to  the  delamination  being  located  at  the 
midsurface  and  the  small  ratio  of  bending  strain  to 
membrane  strain,  snap-through  will  not  occur  on  any  of  the 
panels.  Sample  calculations  of  the  non-dimensionalized 
terms  are  presented  in  Appendix  C. 

STAGSC-1  Theory 

STAGSC-1  (Structural  Analysis  of  General  Shells)  is  an 
energy  based,  finite  element  computer  program.  The  code, 
developed  by  Lockheed  Palo  Alto  Research  Laboratory,  was 
developed  to  analyze  shell  structures  under  various 
loadings  and  boundary  conditions.  The  eigenvalue 
solutions  used  in  this  thesis  were  obtained  from  a  linear 
analysis  utilizing  the  QUAF  411  element.  The  grid  size 
for  each  panel  was  1  inch  by  1  inch.  For  further 
information  on  STAGSC-1  Linear  Analysis  and  grid  modeling, 
the  reader  is  referred  to  references  [5]  and  [18]. 


III.  Manufacturing  and  Experimental  Procedures 


Panel  Manufacturing 

The  panels  used  in  the  compression  tests  were 
constructed  from  Hercules  AS4/3502  Graphite/ Epoxy  12  inch 
tape.  An  initial  set  of  two  flat  panels,  one  [0/90] 2s 
lay-up  and  one  [ 0/-45/ +45/ 90 ] s  lay-up,  were  constructed 
to  determine  a  suitable  material  to  cause  a  de lamination . 
Six  material  combinations  were  tested  in  each  of  the 
panels:  mylar  (t  =  1  mil),  mylar  with  release,  backing 

paper  (t  =  5  mil),  backing  paper  with  release,  non-porous 
teflon  (t  =  3  mil) ,  and  nylon  peel  ply  (t  =  3  mil) .  The 
release  used  was  RAM  225  and  it  provided  a  further 
safeguard  to  prevent  bonding.  After  the  curing  cycles, 
each  panel  was  examined  by  a  C-Scan  to  confirm 
delamination  had  occurred  and  also  to  insure  the 
de laminations  had  not  moved  during  the  lay-up  and  curing 
process.  The  C-Scan  results  indicated  all  six  materials 
had  caused  de laminations  and  none  of  the  materials  shifted 
or  moved  during  the  manufacturing  process.  As  a  further 
check,  one  panel  was  cut  length-wise  through  the 
de laminations  (Fig.  6)  and  notches  were  cut  from  each 
delaminated  section.  Again,  both  sides  of  the  panel  fell 
apart  at  the  delamination  indicating  a  complete  lack  of 
bonding  around  each  material.  The  material  chosen  for  the 
curved  panel  de laminations  was  the  mylar  with  release 


because  it  had  the  smallest  thickness  of  the  four 
materials . 

All  curved  panels  were  hand  laid  on  a  curved  plate 
with  the  required  12  inch  radius  of  curvature.  Mylar 
templates  (Fig.  7)  were  cut  for  each  mold  in  order  to 
place  the  delamination  material  at  the  proper  location 
(Fig.  8).  Following  the  lay-up,  the  panels  were 
transferred  to  the  steel  molds  (Fig.  9).  In  the  molds  the 
panels  were  prepared  for  the  autoclave  and  sealed  in  a 
nylon  vacuum  bag  (Fig.  10).  The  entire  lay-up  was  cured 
according  to  the  following  cure  cycle: 

(1)  Apply  full  vacuum  (2.42  psia  minimum) 
to  sealed  lay-up 

(2)  Heat  to  270°  F  in  45^5  minutes 

(3)  Hold  270°  and  full  vacuum  for  15  minutes 

(4)  Pressurize  autoclave  to  85  psi 

(5)  Hold  at  270r  and  85  psi  for  45  minutes 

(6)  Heat  to  350°  F 

(7)  Hold  at  350° F  and  85  psi  for  15  minutes 

(8)  Vent  lay-up  bag  vacuum 

(9)  Hold  at  350° F  and  85  psi  for  105  minutes 

(10)  Cool  to  below  150° F  in  60  minutes 
maintaining  85  psi 

(11)  When  below  150° F,  vent  pressure 
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Midsurface  Delaminations 
Figure  6 
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Template  to  Position  Delamination 


Figure  7 


Midsurface  Circular  Delaminations 


Figure  8 


After  curing.,  the  panels  were  checked  by  C-Scan 
techniques  to  insure  no  voids  were  present.  In  addition, 
a  hand-held  ultrasonic  test  devise  was  used  on  each  panel 
to  accurately  mark  the  delamination  location,  Fig.  11. 
These  panel  markings  were  permanent  in  order  to  use  as  a 
reference  comparison  after  the  panel  buckling  phase  of  the 
experiment  was  carried  out. 

Each  cured  lay-up  was  cut  into  four  panels  of  the 
desired  experimental  size.  The  fcest  panels  were  cut  on  a 
specially  jigged  radial  arm  saw  with  a  diamond  tipped 
blade.  After  each  panel  was  cut,  a  series  of  8  length,  4 
chord,  and  12  thickness  measurements  were  taken  .  The 
thickness  measurements  showed  the  panel  variation  was  less 
than  1.0  percent  for  any  panel.  The  length  dimensions  are 

very  critical  to  insure  an  even  load  distribution 

0 

throughout  the  test.  From  this  thesis  and  previous 
experiments  conducted  with  the  compression  [5,  9],  it  was 
determined  that  variations  greater  than  .01  inch  on  a  12 
inch  chord  and  .007  inch  on  an  8  inch  chord  were 
unacceptable.  All  chord  measurements  were  acceptable 
since  a  quarter  inch  gap  was  available  on  the  vertical 
side  supports. 
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Hand-held  Ultrasonic  Equipment 
Figure  11 
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Experimental  Set-Up 

The  test  device  used  in  this  thesis  required  a 
modification  of  a  current  compression  testing  system.  This 
device  has  a  30K  lb  capacity  and  was  modified  with  fixtures 
to  hold  the  curved  panels.  Boundary  conditions  assumed  to 
be  provided  by  the  test  set-up  are  a  clamped  top  edge 
(u  =  v  =  w  =  w  x  =  0) ,  clamped  bottom  edge 
(u  =  free,  v  =  w  =  w/X  =  0) ,  and  simple  supports  along 
both  vertical  sides  (u  =  v  =  w;y  =  free,  w  =  0) . 

Both  top  and  bottom  plates  of  the  test  fixture  are 
bolted  directly  to  the  cross  heads  of  the  compression 
machine.  The  load  is  applied  to  the  bottom  cross  head  by  a 
hydraulic  ram.  Clamped  boundary  conditions  are  obtained  by 
a  set  of  blocks  forced  against  the  top  and  bottom  edges 
with  a  series  of  screws.  The  screws  are  tightened  with  a 
screwdriver  to  eliminate  any  panel  movement.  Center  lines 
are  marked  on  both  top  and  bottom  test  fixtures  to  insure 
the  panel  is  properly  centered. 

Simply  supported  side  boundary  conditions  are  provided 
by  two  bars  clamped  in  each  side  support.  The  edges  of  the 
bars  are  thin  knife  like  edges  that  allow  the  panel  to 
rotate  in  the  w^y  direction.  Figure  12  shows  the  side 
supports  with  the  bars  and  plates  used  to  clamp  the  bars  in 
the  side  supports.  When  inserting  the  panel,  a  series  of 
five  set  screws  are  tightened  on  each  side  of  the  support. 
These  set  screws  are  finger  tightened  until  the  knife  edges 
just  contact  the  panel. 
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In  all  tests  the  top  and  bottom  of  the  panel  were 
clamped  first,  then  the  simple  supports  were  installed  and 
the  set  screws  finger  tightened.  A  1/4  inch  space 
occurred  between  the  knife  edges  and  the  bottom  supports 
to  allow  the  panel  to  compress.  By  adjusting  the  knife 
edges,  this  1/4  gap  is  placed  at  the  bottom  to  achieve 
buckling  symmetry  [5].  Fig.  13  shows  the  complete 
experimental  set  up  with  panel  installed. 

When  conducting  buckling  experiments,  two  problems 
have  constantly  plagued  investigators,  those  of  proper 
load  introduction  [19],  and  the  effects  of  boundary 
conditions  changing  during  the  experiment  [14].  Both  of 
these  problems  were  encountered  during  the  experiments  for 
this  thesis.  However,  by  carefully  following  the  above 
procedure  for  panel  installation,  the  effects  of  these  two 
problems  were  minimized. 

Instrumentation  consisted  of  four  axial  strain  gages 
and  an  LVDT  (Linear  Variable  Differential  Transducer). 

The  LVDT,  visible  on  the  right  side  of  the  test  fixture  in 
Fig.  14,  was  mounted  on  the  cross  heads  to  measure 
end-shortening  of  the  panel.  Two  of  the  strain  gages  were 
mounted  two  inches  from  the  bottom  of  the  panel  and  three 
inches  from  the  panel  center  line  (two  inches  from  panel 
omn-r  line  on  panels  with  an  eight  inch  circumference)  to 
insure  uniform  loading.  Two  additional  gages  were 
mounted  back-*- o-back  at  the  geometric  center  of  the 
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confirmed  by  examining  the  experimental  results.  The 
reader  will  recall  the  two  strain  gages  placed  equidistant 
from  the  center  line  to  determine  whether  the  load  was 


applied  symmetrically.  An  unsymmetric  load  application 
resulted  in  one  gage  showing  more  compressive  strain  than 
the  other  gage.  Since  the  side  showing  more  strain  also 
was  taking  more  of  the  load,  the  first  buckle  occurred  on 
this  side  of  the  panel.  By  further  increasing  the  load, 
the  other  side  of  the  panel  began  taking  the  load  and 
would  then  buckle.  This  phenomena  was  confirmed  by 
comparing  buckl ing/ pos tbuckl ing ,  strain  gage  output,  and 
buckling  patterns.  However,  analyzing  the  varying 
boundary  conditions  and  the  postbuckling  problem  are 
beyond  the  scope  of  this  thesis. 

Plots  of  load  versus  end-shortening  for  the  analytical 
and  experimental  values  are  shown  in  Figures  16-19  through 
the  first  buckle  only.  These  four  figures  indicate  the 
reduction  in  stiffness  for  each  panel  as  the  delamination 
size  increases.  Also,  when  comparing  panels  with  equal 
aspect  ratios  and  different  ply  orientations,  it  is  shown 
the  [ 0/ -45 / + 45/ 90 ] s  panels  have  a  greater  reduction  in 
stiffness  as  delamination  size  increases  when  compared  to 
the  (0/90l2s  panels. 
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Panel  ID 

Buckling  Load  (lb) 

End-Shortening  (in) 

Q201 

5574.6 

.0842 

Q221 

5027.4 

.0833 

Q222 

4663.. 3 

.0801 

Q242 

2527.1 

.0493 

C201 

4275.0 

.0732 

C221 

4181.9 

.0719 

C242 

3453 . 9 

.0678 

Q802 

6246.9 

.0700 

Q822 

4135.6 

.0586 

Q841 

4157.0 

.0529 

Q842 

3796.3 

.0555 

C802 

5607.7 

.0528 

C821 

4825.5 

.0480 

C822 

4610.0 

.0456 

C841 

4038.0 

.0398 

C84  2 

4485 . 2 

.0424 

Experimental  Results 
Table  I 
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providing  the  desired  boundary  conditions.  The  frequent 
unsymmetric  buckling  patterns  and  the  postbuckling 
strength  exhibited  by  the  panels  led  to  the  assumption 
that  the  vertical  side  supports  were  providing  more 
restraints  than  the  desired  simple  support.  To  verify 
this  assumption,  STAGSC-1  modeling  was  accomplished  by 
varying  the  vertical  side  boundary  conditions.  The 
results  of  numerous  previous  buckling  experiments  have 
confirmed  that  an  analytical  buckling  load  with  no  panel 
imperfections  is  higher  than  an  experimentally  determined 
load  for  the  same  panel.  Therefore,  the  initial  modeling 
consisted  of  simply  supported  {u  =  v  =  w(y  =  free,  w  = 

0)  boundary  conditions  and  then  eliminating  certain 
degrees  of  freedom.  From  the  simple  support  conditions, 
first  u  was  fixed  and  then  both  u  and  v  fixed.  To  achieve 
an  analytical  bifurcation  load  higher  than  the 
experimental  value,  it  was  necessary  to  have  both  u  and  v 
fixed  along  the  vertical  sides.  These  are  the  analytical 
values  presented  in  Table  II.  In  addition  to  increasing 
the  initial  buckling  load,  a  stiffer  panel  will  exhibit- 
greater  postbuckling  strength  as  was  evidenced  in  these 
experimental  tests.  To  properly  analyze  this  phenomena, 
one  would  have  to  vary  the  boundary  conditions  as  the  load 
is  applied  and  the  panel  becomes  progressively 
restrained.  To  further  complicate  the  problem,  the 
restraints  are  not  applied  symmetrically.  This  finding  is 
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Table  I  shows  the  experimental  results  for  the 
panels.  The  reader  will  notice  that  of  the  24 
manufactured  panels,  only  16  are  included  in  the  results. 
Of  the  remaining  eight  panels,  seven  were  rejected  due  to 
exceeding  the  length  dimension  tolerance.  These 
tolerances  of  .01  inch  variation  on  panels  with  AR  =  1  and 
.007  inch  variation  on  panels  with  AR  =  1  1/2  were 
reconfirmed  after  examining  the  experimental  output. 
Exceeding  these  tolerances  prohibited  an  even  load 
distribution  through  the  panel,  verified  by  the 
unacceptable  differences  between  the  two  strain  gages 
mounted  to  monitor  load  introduction  and  what  appeared  to 
be  shear-induced  buckling  patterns. 

Analytical /Experimental  Comparison 

To  this  author's  knowledge,  there  are  no  analytical 
expressions  or  computer  coded  models  that  will  solve  the 
problem  of  a  curved  composite  panel  with  a  delamination 
located  at  the  midsurface.  Therefore,  in  order  to 
establish  limits  on  the  experimental  results,  the  STAGSC-1 
code  was  used  to  model  the  panels  with  no  delamination  and 
establish  an  upper  bound  on  the  buckling  loads.  These 
analytical  results  and  final  experimental  values  are 
presented  in  Table  II. 

It  became  apparent  during  the  initial  stages  of  the 
compression  testing  that  the  test  fixture  was  not 
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versus  panel  end-shortening  and  load  versus  back-to-back 
strain  gages  are  used  to  determine  the  panels'  buckling 
load.  The  third  plot  was  load  versus  the  two  remaining 
strain  gages  and  is  used  to  insure  the  load  is  evenly 
applied  to  the  panel. 

The  use  of  a  load  versus  deflection  (or  end¬ 
shortening)  plot  to  determine  buckling  load  was  explained 
in  reference  to  Fig.  2.  Another  technique  to  determine 
buckling  loads  is  known  as  a  strain  reversal  technique. 
This  method  is  shown  in  Fig.  15  and  requires  the  output  of 
back-to-back  strain  gages.  As  the  load  is  initially 
applied  and  begins  to  increase,  both  strain  gages  register 
the  same  amount  of  compression.  The  gages  continue  to 
read  the  same  increase  in  strain  until  the  panel  buckles. 
At  this  point,  the  strains  of  the  respective  gages 
increase  dramatically  in  both  tension  and  compression 
indicating  the  physical  presence  of  a  buckle.  One 
limitation  to  this  technique  is  that  it  requires  a  prior 
knowledge  of  the  panel's  buckling  pattern  to  insure  the 
back-to-back  gages  are  located  at  the  initial  buckle.  The 
back-to-back  gages  on  these  panels  were  located  at  the 
geometric  center  of  the  panel,  which  is  also  the  center  of 
the  de laminations .  Since  the  initial  buckle  always 
occurred  at  the  de lamination ,  the  strain  reversal 
technique  yielded  the  same  buckling  load  as  the  load 
versus  deflection  plot. 


IV.  Results  and  Discussions 


Panel  Identification 

The  experimental  panels  are  identified  using  the 
following  system: 

Panel  1  or  duplicate  Panel  2 

Delamination  Size 
0  =  no  delamination 
2=2  inch  diameter  delamination 
4=4  inch  diameter  delamination 

Aspect  Ratio 

2  =  12  inch  length  X  12  inch  chord 
(AR  =  1) 

8  =  12  inch  length  X  8  inch  chord 
(AR  =  1  1/2) 

Ply  Orientation 

Q  =  [0/-45/+45/90] s ,  quasi-isotropic 
C  =  [0/90] 2s  •  cross-ply 

A  three  digit  code  is  used  to  identify  the  average  results 

of  both  duplicate  panels.  As  an  example,  Q22  means  the 

average  values  for  panels  Q221  and  Q222. 

Data  Output 

As  stated  previously,  all  data  was  stored  on  tape  and 
printed  out  at  a  later  time.  The  data  from  all 
instrumentation  was  referenced  to  a  computer  reference 
time  and  hence  the  output  from  any  instrument  cb~  el 
could  be  compared.  In  addition  to  the  printed  data,  three 
computer  generated  plots  were  available.  Two  plots,  load 


panel.  This  instrumentation  provided  the  data  to 
determine  buckling  loads  by  two  techniques: 

1 .  load  versus  end-shortening 

2.  strain  reversal. 

Test  Procedure 

Once  the  panel  was  installed,  the  LVDT  and  all  strain 
gages  were  electronically  zeroed.  The  load  was  applied  at 
the  rate  of  .055  inches  per  minute  cross  head  speed  and 
data  was  recorded  ten  times  per  second  on  a  VAX  11/780 
computer.  Loading  was  held  constant  through  the  initial 
buckle  and  beyond  in  order  to  get  a  portion  of  the 
postbuckling  curve.  The  buckling  pattern  was  highlighted 
with  a  silver  pencil  and  the  pattern  photographed.  All 
data  was  stored  in  a  real  time  mode  and  printed  out 
later.  Data  from  the  load  cell  and  the  four  strain  gages 
was  also  displayed  on  a  CRT  during  the  tests  to  monitor 
load  and  load  introduction  to  the  panel. 


Complete  Experimental  Set-Up 
Figure  14 
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Buckling  Patterns 

When  installing  the  panel  in  the  test  fixture,  it  was 
found  that  by  leaving  the  1/4  inch  of  free  space  at  the 
loading  edge  of  the  panel  the  buckling  pattern  would  be 
symmetric  [5]-  The  gap  at  the  end  where  the  load  was 
applied  allowed  for  a  more  even  stress  distribution.  This 
symmetric  buckling  pattern  can  be  seen  in  Fig.  20,  a  panel 
with  no  delamination.  But,  the  buckling  pattern  differed 
when  compared  to  panels  that  had  de laminations .  All 
panels  with  de laminations  exhibited  the  initial  buckle  at 
the  delamination  location.  However,  the  buckling  patterns 
differed  after  the  initial  buckle  appeared  and  the  load 
was  increased.  Figures  21  and  22  represent  two  different 
post-buckling  patterns.  Fig.  21  is  representative  of  a 
symmetric  pattern  where  the  initial  buckle  grew  about  the 
delamination  and  the  pattern  remained  symmetric  as  two 
buckles  appeared  in  the  lower  corners.  In  contrast, 

Fig.  22  shows  a  buckling  pattern  that  grew  from  the 
initial  buckle  to  a  shear-type  buckling  pattern  indicating 
a  variable  change  in  the  boundary  conditions  on  the 
vertical  supports. 


Figure  20 


Q841  Buckling  Pattern 
Figure  21 
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Q822  Buckling  Pattern 
Figure  22 
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Postbuckling  Panel  Inspection 


The  reader  will  recall  the  purpose  of  this  thesis  is 
to  experimentally  determine  the  buckling  load  of  a 
delaminated  panel  where  the  delaminated  layer  maintains 
full  contact  with  the  panel-  Visual  observation  of  the 
buckling  pattern  for  each  test  confirmed  snap-through  did 
not  occur  on  any  panel.  To  insure  no  delamination  growth 
or  any  other  panel  damage  occurred  during  the  buckling, 
two  techniques  were  employed-  First,  each  panel  was 
scanned  using  the  same  hand-held  ultrasonic  transducer 
(Fig.  11)  used  to  locate  the  delamination  after  the  panels 
were  manufactured.  Since  each  delamination  had  been 
permanently  marked,  it  was  a  simple  procedure  to  check  for 
any  delamination  growth.  By  moving  the  transducer  across 
the  delamination,  the  edges  were  located  and  found  to 
exactly  coincide  with  the  previous  delamination  markings 
for  all  panels.  Thus,  it  was  verified,  no  delamination 
growth  occurred  for  any  panels. 

The  second  method  involved  the  use  of  x-rays  to  check 
for  any  matrix  cracking.  This  technique  required  the 
introduction  of  a  penetrant,  tetrabromoethane ,  into  the 
delaminated  area  and  insuring  the  penetrant  completely 
covered  this  damaged  area.  One  panel  from  each  duplicate 
set  was  used  for  an  x-ray  by  drilling  a  1/2  inch  hole 
through  the  de lamination .  The  penetrant  was  then  brushed 
around  the  hole  and  allowed  to  dissipate  throughout  the 


delaminated  area.  An  x-ray  will  show  the  areas  covered  by 
the  penetrant  as  darker  areas  contrasted  with  the  normal 
panel.  Fig.  23  is  a  typical  x-ray  pattern.  One  will 
notice  the  different  shaded  areas  in  the  delamination 
caused  by  three  applications  of  the  penetrant.  However, 
of  more  importance  is  the  edges  of  the  delaminated  area. 
The  penetrant  did  not  seep  past  the  well  defined  circular 
delamination  indicating  no  panel  damage,  matrix  cracking, 
or  delamination  growth  occurred. 


Ply  Or  Mentation/ Aspect  Ratio  Effects 


V 


Table  II  shows  all  analytical  and  experimental 
buckling  loads,  as  well  as,  percent  decrease  in  the  loads 
for  each  size  delamination.  As  was  expected,  the  larger 
the  size  of  the  delamination,  the  less  the  buckling  load. 
The  [ 0 /  —  4 5 / +45/ 90] s  orientation  had  the  highest  buckling 
value  for  panels  of  equal  aspect  ratio  with  no 
delamination,  and  also  showed  the  greatest  decrease  in 
experimental  values  for  each  size  de lamination .  This 
result  was  also  expected  due  to  the  lay-up  being  more 
dependent  on  the  _+  45  orientations  in  both  the  extensional 
and  bending  stiffness  matrices.  Values  for  all  stiffness 
matrices  are  presented  in  Appendix  B.  Figures  24  and  25 
graphically  show  the  higher  experimental  buckling  values 
of  the  quasi-isotrooic  lay-up  with  no  de lamination ,  and 
the  greater  decrease  in  the  experimental  buckling  load  of 
this  lay-up  as  delamination  size  is  increased,  when 
compared  to  the  cross-ply  orientation.  The  vertical  axis 
representing  the  buckling  loads  is  normalized  by  dividing 
each  experimental  load  by  the  largest  value,  which  is  the 
buckling  load  of  the  [ 0/ -4 5/ + 45/ 90 ] s  orientation  with  no 
de laminat ion .  The  decrease  in  the  load-handling 
capability  of  each  panel  as  delamination  size  increases 
corresponds  to  the  knockdown  factors  given  in  Table  II. 
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Non-Dimens iona 1 i zed  Experimental  Values 
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Figure  24 
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Non-Dimens lonalized  Experimental  Values 
versus  De lamination  Diameter  AR  =  1  1/2 


The  reader  will  note  from  Figures  24  and  25,  the 
^  [0/90] 2s  panels  become  stiffer  at  a  delamination 

diameter  of  approximately  3  inches  ori  the  panels  with  AR  = 
1  (Fig.  24),  and  a  delamination  diameter  of  approximately 
1  inch  on  the  panels  with  AR  =  1  1/2  (Fig.  25).  At  these 
points,  the  buckling  problem  becomes  one  of  the  axial 
stiffnesses  restraining  the  panel,  rather  than  the  bending 
I  stiffnesses  that  restrain  a  panel  with  no  or  small 

de laminations .  Figures  16  and  17  confirm  the  [0/90] s 
panel  becomes  stiffer  at  a  3  inch  delamination  for  panels 
.  with  an  AR  =  1.  By  choosing  an  Nx  and  interpolating  to 

a  3  inch  de lamination ,  one  will  find  the  [0/90] 2s  panel 
is  stiffer  since  it  allows  less  end-shortening.  Comparing 
f  Figures  18  and  19  and  interpolating  to  a  1  inch 

delamination  confirms  the  [0/90] 2s  becomes  stiffer  at 
the  1  inch  or  larger  delamination  for  an  AR  =  1  1/2. 

« 

De lamination/Cutout  Comparison 

Previous  work  by  Janisee  [9]  determined  the  analytical 
collapse  loads  for  curved  panels  with  square  cutouts.  His 
results  of  a  [0/-45/ +45/ 90] s  orientation  and  an  AR  =  1 
are  compared  to  the  same  panel  with  de laminations .  The 
analytical  results  are;  Nx  =  252.1  Ib/in  for  a  2X2  inch 
cutout  and  Nx  =  113.2  for  a  4X4  inch  cutout  for  a  panel 
with  simple  support  conditions  along  the  vertical  sides. 
These  values  result  in  reduction  of  56.1%  and  80.3%  when 
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compared  to  the  analytical  bifurcation  load  with  no 
defects,  and  represent  a  definite  lower  bound  to  the 
delaminated  panel  since  the  simple  support  boundary 
conditions  yield  a  panel  that  is  less  stiff  then  boundary 
conditions  provided  by  the  experimental  fixture.  The  2X2 
inch  cutout  load  is  37.6%  less  than  the  load  with  a  2  inch 
de lamination ,  while  the  4X4  inch  cutout  load  is  46.2%  less 
than  the  4  inch  delamination  load. 

As  was  previously  discussed  in  the  Analytical/Experi¬ 
mental  Comparison  section,  a  variety  of  boundary  condi¬ 
tions  along  the  vertical  sides  were  modeled  in  the 
STAGSC-1  code.  These  results  confirmed  the  experimental 
fixture  did  not  provide  the  desired  simple  support  condi¬ 
tions  along  the  vertical  sides,  and  also  established  an 
upper  bound  to  the  problem  of  determining  the  buckling 
load  of  a  delaminated  panel.  By  restraining  degrees  of 
freedom  along  the  vertical  edges,  a  stiffer  panel  with  a 
higher  bifurcation  load  resulted. 

Three  different  sets  of  boundary  conditions  along  the 
vertical  sides  were  considered  with  the  STAGSC-1  computer 
code.  First,  simple  support  conditions  were  modeled 
allowing  u,v  and  rotation  about  the  x-axis  to  be  free 
along  the  entire  vertical  sides.  A  second  choice 
restrained  u  and  v  movement,  allowing  only  rotation.  The 
third  set  of  boundary  conditions  allowed  rotation  along 
the  entire  sides  and  u  and  v  free  along  the  top  one  inch 


of  the  panel  where  the  load  is  applied.  The  purpose  of 
this  third  choice  was  to  model  the  changing  boundary 
conditions  of  the  experimental  fixture  as  the  load  is 
applied  and  the  set  screws  tighten  the  knife  edges, 
restraining  panel  movement  (see  Experimental  Set-up).  The 
results  of  these  three  sets  of  boundary  conditions  are: 


U  ,  V 

fixed 

Nx 

=  574.0 

lb/in 

U  ,  V 

free  along 

top 

one  inch 

NX 

=  538.1 

lb/  in 

U  ,  V 

free 

NX 

=  507.9 

lb/  in 

The 

results  of 

this 

analysis 

are 

plotted 

in  Fig 

with  the  vertical  axis,  buckling  load,  normalized  to  the 
analytical  buckling  load  of  a  perfect  panel  with  u  and  v 
fixed  along  the  vertical  sides,  Nx  =  574.0  lb/in.  The 
three  values  corresponding  to  1.0,  .94  and  .88  represent 

the  loads  of  the  varying  boundary  conditions;  u,v  fixed, 
u,v  free  along  top  one  inch,  and  u,v  free,  respectively. 
Each  of  these  loads  is  for  a  panel  with  no  de laminations 
or  cutouts,  hence  they  provide  an  upper  bound  to  the 
delaminated  panel  with  various  boundary  conditions. 


l.o 


.8 


0  2 

Delamination  Diaraeter/Cutout  Size 
( inch) 

X  -  u#v  fixed 

■4-  -  urv  free  along  top  one  inch 

S-  u#v  free 

-  delamination 
Q  -  cutout 
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Non-Dimensionalized  Buckling  Load  versus  Defect  Size 

[ 0 /  —  4  5 /  +  45/90] s  AR  =  1 

Figure  26 
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V.  Conclusions 


Based  on  the  analytical  and  experimental  results 
presented  in  this  thesis,  the  following  conclusions  can  be 
made : 

1.  To  insure  valid  experimental  results,  the 
researcher  must  carefully  choose  the  proper  test 
device.  When  conducting  experiments  on  axially 
compressed,  curved  composite  panels,  the  test 
fixture  must  insure  an  even  load  application  to 
the  panel  and  provide  the  desired  boundary 
conditions . 

2.  The  STAGSC-1  linear  analysis  computes  a 
bifurcation  load  for  a  perfect  (no  eccentricities 
or  defects)  panel.  This  load  provides  an  upper 
bound  for  the  experimental  buckling  load  of  a 
panel  with  the  same  ply  orientation,  aspect 
ratio,  and  boundary  conditions.  A  lower  bound  to 
the  experimental  buckling  load  of  a  panel  with  a 
defect  is  provided  by  the  STAGSC-1  analytical 
results  of  a  panel  with  a  cutout.  These  two 
results  bound  the  experimental  value  of  a 
delaminated  panel  with  a  random  size  and  shape 
defect . 
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3.  The  process  used  to  manufacture  the  experimental 
panels  will  effect  the  experimental  results. 
Panel  variations  in  the  length,  width,  and 
thickness  dimensions  need  to  be  accounted  for 
when  modeling  the  panel  for  an  analytical 
analysis . 

4.  Putting  a  delamination  at  the  midsurface  of  a 
curved  panel  allows  for  testing  of  a  specimen 
where  the  delaminated  layer  does  not  separate 
from  the  panel  and  snap-through  does  not  occur. 

5.  The  delamination  size  effects  the  load-handling 
capability.  The  larger  the  delamination,  the 
lower  the  experimental  buckling  load  of  the 
pane  1 . 

6.  The  quas i- isotropic ,  [0/-45/  +  45/90] s  , 

orientation  exhibited  a  higher  analytical  and 
experimental  buckling  load  than  the  cross-ply, 
[0/90]2s»  orientation  for  panels  of  the  same 
aspect  ratio  and  no  de lamination .  However,  for 
panels  with  the  same  size  de lamination ,  the 
quas i- iso  tropic  panels  had  the  greatest  decrease 
in  buckling  loads  and  for  certain  delammation 
sizes  the  quasi-isotropic  yielded  lower  loads 
than  the  cross-ply  panels. 
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Similar  calculations  for  each  experimental  panel 


showed  no  value  of  a  to  exceed  3-347.  Therefore, 
analytical ly ,  no  panel  would  experience  snap-through  or 
delamination  growth.  This  result  was  experimentally 
verified . 


Since  the  experiments  for  this  thesis  applied  a  pur® 
compressive  load  with  bending  strains  only  induced  by  the 


panel  curvature,  a  comparison  was  made  to  the  minimum 
value  of  p computed  by  Fei  and  Yin  which  was  .01.  By 
computing  the  value  of  a  for  each  experimental  panel,  it 
was  confirmed  snap-through  would  not  occur  so  long  as  the 
computed  value  for  a  did  not  exceed  3.347.  Equation  20 
is  the  expression  for  a.  The  value  for  uwas  the 
effective  laminate  property  computed  in  Appendix  B.  The 
strain  in  the  delaminated  layer,  t,  is  the  strain 
indicated  by  the  back-to-back  strain  gages  centered  on 
the  delamination.  The  delamination  depth  ratio,  a/h,  is 
determined  from  each  size  delamination  and  a  delamination 
depth  of  .022  inch. 

Panel  Q822 

v  =  .307,  (=  329.191  X  10~6 ,  a  =  1  inch,  h  =  .022  inch 
a  =  [12 ( 1-U2  ) e ] 1/2  a/h  = 

[12  (1-.  3072)  (329.191  x  i0-6 ) ]  1/2  (1/ . 022) 
a  =  2.719 

.  J 

j 

Panel  C242 

'  ^ 

j 

v  -  .04,  «=  84.319  X  10-6,  a  =2  inch,  h  =  .022  inch  j 

a  =  [12 ( 1- .  0  4  2 )  (84.319  X  10~6 ) ]  2  (  2/ .  022 )  J 

a  =  2.889 


(IOC) 


a  a  a 

i«  +  l/2  J  M2r2dr  =  -  Jt rdr  +  1/2  y~(w')2di 


Applying  ( 6  C )  to  (7C)  yields: 


E6r  =  <Tr-  ucr^  =  |-P+  u(rP)  ' 


/h  ( 11C ) 


After  integration  by  parts  twice  on  the  right  hand  side  of  (IOC) 
and  substituting  ( 1 1C ) : 

(1-u)  P  (a)  +  a  P'  (a)  =  Eh(*  +  M2a2/6)  (12C) 

The  system  of  four  equations  (7C) ,  (8C) ,  (9C) ,  (12C)  completely 

determine  the  functions  <*>(r)  and  P(r). 

However,  by  non-dimensionalizing  the  parameters,  the 
solution  to  the  complete  buckl ing/postbuckl ing  problem  is 
reduced  to  the  three  factors:  p,  a,  Ga^.  The 
phenomena  of  snap-through  is  described  by  p and  a  only, 
since  their  analysis  found  that  so  long  as  the 
delaminated  layer  maintains  full  contact  with  the  main 
body  of  the  plate,  delamination  growth  can  not  preceed 
snap- through •  Therefore  the  results  of  Fei  and  Yin  can 
be  used  to  determine  which  of  the  three  states  of  the 
buckling/postbuckl ing  problem  a  panel  will  experience. 

Their  numerical  results  indicated  that  a  minimum  value 
of  pis  .01, which  corresponds  to  a  minimum  a  of  3.347. 
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Let  P{r)  be  the  compressive  radial  force  per  unit 
circumf erential  arc  length,  and  w(r)  the  transverse 
deflection.  With  the  applied  force,  the  delaminated  layer 
will  either  maintain  contact  with  the  main  body  of  the  plate 
or  buckle  away  from  the  plate,  defined  as  snap-through.  In 
the  first  case,  the  deformed  shape  of  the  plate  is: 

1/fX  =  1  /fxo  -  H/2  ;  H  =  t-h  (5C) 

The  condition  of  snap-through  is  represented  by  <t>: 

<t>  (r)  =  w'  (r)  /r  (60 

The  P  and  <t>  satisfy  a  system  of  two  second  order  non-linear 
differential  equations: 

r“3  (r3P'  )  '  =  Eh<o2/2  (7C) 

r-3 (r3o' ) ’  =  P0/D  (8C) 

for  r£  a  and  the  boundary  conditions  at  the  center  of  the 
layer: 

P'  (0)  =  0,  0'  (0)  =  0  (9C) 

Countinuity  of  slope  along  the  boundary  of  the 
delamination  requires  w'(a)  =  v'(a).  Also,  by  comparing  the 

meridional  arc  lengths  of  the  unbuckled  and  buckled 
delamination  it  is  found: 


r 


Appendix  C 

» 

Analytical  Non-Dimensionalized  Values 

The  results  of  the  Fei  and  Yin  [13]  were  used  to 
S'  confirm  the  delaminated  layer  would  not  snap-through . 

Their  results  presented  snap-through  states,  which  give 
the  required  delamination  radius,  a,  to  cause  snap-through 
®  for  selected  values  of  p,  ratio  of  bending  deformation  to 

membrane  deformation.  The  analysis  uses  Bessel  functions 
to  express  the  slope  and  curvature  of  the  deformed  middle 
1  surface  in  terms  of  the  radial  coordinate,  r,  and  the 

ratio  of  membrane  force  to  the  plate  bending  stiffness,  a. 

•  v’  (r)  =  A  J1  (  x  r)  (1C) 

v"  (r)  =  A  A  [J0A  (  A  r)-(l/Ar)J1(  A  r)  ]  (2C) 

*  By  applying  the  power  series  expansions  of  Bessel  functions, 
and  assuming  (  A  a«l)  ,  a  small  membrane  force  in  the  plate 
the  equations  simplify  to: 


V  (a)  =  p0  a  (3C) 

v"  (a)  =  Mo  (40 


where  u  =  aX/2  is  the  curvature  of  the  middle  surface  at 
the  center  of  the  plate. 
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[ 0/ -45/ +45/ 90 ]  g; 


.361 

•  111 

0 

.111 

.361 

0  106 

0 

0 

.125 

Bij  =  0 


79.867 

12.502 

-4.884 

12.502 

21.162 

-4.884 

-4.884 

-4.884 

14.392 

Ex  =  Ey  =  8.172  X  106  psi 

°xy  ~  uyx  =  .307 

3xy  =  3.125  X  106  psi 


[0/90] 

Aij  = 


Bij  « 


.453  .018 
.018  .453 
0  0 


0 

0 

.033 


106 


0 


Bij 


80.117  2.470 

2.470  40.979 

0  0 


0 
0 

4.3  561 


,10 


-x  =  Ey  =  11.307  X  106  psi 
l>Xy  -  ~Uyx  =  *0  40 

=  .825  X  106  psi 


Appendix  B 


Extensional  Stiffness 


’  A11 

a12 

a16 

Aij  - 

a12 

a22 

a26 

.  a16 

a26 

a66  . 

Coupling  Stiffness 

B11 

b12 

b16 

Bij  = 

b12 

b22 

b26 

.  B16 

b26 

b66  , 

Note:  For  a  symmetric 

lay-up 

<  Bi j 

=  0 

Bending  Stiffness 

Dll 

d12 

Die 

Dij  = 

d12 

d22 

d26 

.  D16 

d26 

d66 

Effective  Laminate  Properties  [16] 

Longitudinal  Youngs'  Modulus:  Ex  =  (Ai]_A22~A12^ 
Transverse  Youngs'  Modulus:  Ey  =  (  A;liA22-a122 

Longitudinal  Poissons'  Ratio:  uXy  =  ^12/^22 
Transverse  Poissons'  Ratio:  i/yX  =  A^/An 

Shear  Modulus:  GvV  =  Aec/t 


)  /  tA22 
)/tAn 


analysis  is  presented  below. 


Density 

Resin 

|  Weight  (%) 

Volatile 
Weight  (%) 

Void 

Weight  (%) 

Good 

Samples 

1.6211 

'  25.6905 

.0909 

.2366 

Voidy 

Samples 

1.5762 

26.4328 

.2400 

2.8678 

As  can  be  seen  from  the  comparison,  the  voidy  samples  were 
less  dense  in  the  matrix  due  to  the  higher  percentage  of 
voids.  The  material  properties  obtained  from  the  voidy 
samples  are: 


El 

=  20.82  X 

106 

psi 

E2 

=  1.41  X 

106 

psi 

g12 

=  .712  X  106  psi 

y12 

=  .293 

All  the  properties  obtained  from  the  voidy  samples  were 
less  than  those  from  the  good  samples.  The  percent 
decrease  in  values  are:  -  .72%,  E2  -  7.84%,  G^2 

-  12.64%,  v-^2  ~  2.33%.  The  values  of  E^ ,  E2r  and 
i>i2  are  basically  a  measure  of  the  fibers  in  the 
composite  and  do  not  show  a  large  decrease  in  any  of  the 
values.  However,  G^2  is  determined  from  shear  tests  and 
a  composite  transfers  shear  stresses  through  the  matrix. 
The  large  decrease  in  G12  is  then  accounted  for  by  the 
voids  in  the  matrix  unable  to  transfer  the  shear. 
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This  appendix  presents  the  results  of  experimentally 
determined  material  properties  for  two  identical  composite 
lay-ups  cured  under  different  autoclave  cycles.  All 
material  properties  were  calculated  in  accordance  with 
ASTM  standards.  ASTM  D  3039  Test  for  Tensile  Properties 
of  Fiber-Resin  Composites  was  used  to  determine  E]_, 

E2  /  and  l>i2*  ASTM  D  3518  Inplane  Shear  Stress-Strain 
Response  of  Unidirectional  Reinforced  Plastics  was  used  to 
determine  the  shear  modulus  Gi2-  A  total  of  ten  samples 


were  used 

for  each 

test  and  the 

resultant  material 

property 

is  an  average  value  of 

the  ten  samples. 

Ei  = 

20.97  X  106 

psi 

II 

CN 

w 

1.53  X  106 

psi 

Gl2  = 

.815  X  10 

6  psi 

u12  = 

.300 

In  addition  to  the  samples  used  to  determine  the  above 
material  properties,  a  second  set  of  tests  were 
accomplished.  The  samples  for  the  second  test  were  cut 
from  panels  that  were  not  correctly  cured  during  the 
autoclave  cycle.  C-Scan  results  from  these  panels 
revealed  numerous  voids  throughout  the  panels  and  the 
analysis  confirmed  the  improper  resin  and  void  weight.  A 
comparison  of  the  average  values  from  the  chemical 
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Restraining  boundary  conditions  along  the 
vertical  sides  introduce  a  higher  panel  stiffness 
which  resulted  in  larger  buckling  loads. 

The  strain-reversal  and  end-shortening  techniques 
gave  identical  buckling  loads  on  all  panels  with 
a  delamination.  This  is  due  to  the  back-to-back 
strain  gages  located  on  the  delamination  which 
was  the  location  of  the  initial  buckle. 

The  knockdown  factors  for  all  panels  show  the 
same  trends  for  each  ply  orientation,  aspect 
ratio,  and  delamination  size. 
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